Spin-gap phase in nearly-half-filled one-dimensional conductors coupled with phonons 
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Asymptotic properties of nearly-half-filled one-dimensional conductors coupled with phonons are 
studied through a renormalization group method. Due to spin-charge coupling via electron-phonon 
interaction, the spin correlation varies with filling as well as the charge correlation. Depending 
on the relation between cut-off energy scales of the Umklapp process and of the electron-phonon 
interaction, various phases appear. We found a metallic phase with a spin gap and a dominant 
charge-density-wave correlation near half filling between a gapless density-wave phase (like in the 
. doped repulsive Hubbard model) and a superconductor phase with a spin gap. The spin gap is 

produced by phonon-assisted backward scatterings which are interfered with the Umklapp process 
constructively or destructively depending on the character of electron-phonon coupling. 
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I. INTRODUCTION 



Metallic states of doped copper oxides and related materials are extensively studied to clarify the characteristics 
of strongly correlated electron systems. The "spin gap" feature of the underdoped compounds M is one of the most 
peculiar features of them, but its origin is not yet clear. In a model with strong on-site attraction, there would 
i^i - be a temperature range in which local singlet pairs are formed with finite binding energy but the system is still not 
superconducting. In the case of the strong on-site repulsion, however, the spin gap formation must be the consequence 
of a subtle balance. As a fundamental property of strongly correlated materials, it is worth investigating the spin 
excitation gap in metals from various points of view. It is pointed out that a spin gap may play an important role both 
in determining the type of a metal-insulator transition Q and in the mechanism of the high-T c superconductivity. Q| 
It has been pointed out that electron-phonon interaction can drastically change local quasiparticle characters of 
nearly-half-filled electron systems when the on-site repulsion is strong. [|]-[3) In this paper, we investigate asymptotic 
behaviors of correlation functions by integrating out short-distance degrees of freedom in space and time. Repulsive 
interaction and retarded attraction mediated by phonons are treated on an equal footing. Our studies are limited here 
0^ ■ to one dimension where a bosonization technique |t],|8| can be used to derive renormalization group equations and to 
take the filling factor explicitly into account. Attention is paid to tendency for the formation of a spin gap near half 
filling and connections to the insulator phase at half filling and to a possible superconductor phase away from half 
filling. 

Without electron-phonon interaction, spin and charge excitations are separated through bosonization which only 
describes low-energy excitations. Phonons couple with both spin and charge so that they are interfered with each 
other until electron-phonon interaction is integrated out at energy of the order of phonon frequency. This spin-charge 
O ■ coupling at finite energy leads the spin and charge correlations to depend on each other. Technically, it leads the 
equivalent two-dimensional classical system to have interacting Burgers vectors. || 

At half filling, the Umklapp process tends to open a charge gap. Away from half filling, as the energy scale is made 
smaller, the Umklapp process becomes ineffective at finite energy which vanishes in the limit of half filling. Therefore 
the charge excitation is always gapless except half filling. Nevertheless the Umklapp process at finite energy modifies 
the charge correlation and also the spin correlation through the electron-phonon coupling. Thus the spin excitation 
is controlled by the relation between the cut-off energy scale of the Umklapp process and that of electron-phonon 
interaction. In terms of two-dimensional classical systems, the interactions between Burgers vectors originating from 
these two processes have different shapes and different length scales so that they are cut off in different manners. We 
numerically solve renormalization group equations for the correlation exponents and scattering parameters (fugacities 
in the classical system) and find a new phase which is made possible by both of these two processes. 

In this paper, we study a bosonized Hamiltonian in which the standard electronic part is coupled with dispersionless 
phonons so that back and forward scatterings occur exchanging phonons as well as directly between electrons. In the 
electronic part, we consider on-site repulsion. The scattering parameters of electron-electron interaction are taken 
from the lowest-order perturbation with respect to on-site repulsion. This is because the exactly known, fixed-point 
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values of the pure Hubbard model |10[ are no longer valid due to modified scaling properties. For a phonon-assisted 
backward scattering, we take a continuum version of the on-site Holstein coupling |llf| or that of the intersite Su- 
Schrieffer-Heeger (SSH) coupling. |l2| Furthermore, we add a phonon-assisted forward scattering which appears in the 
continuum version of the Holstein model in order to make a singlet superconductor phase possible. It is numerically 
confirmed that this addition to the SSH case does not affect the spin correlation so much. 

The renormalization group equations we obtain deal with the lowest-order perturbations, so that the fixed point of 
either correlation exponent at zero should correspond to the Luther-Emery state [ fL3| with a gap in the corresponding 
channel. At half filling, our results for the critical coupling strengths for the Holstein model and the SSH model 
(without a phonon-assisted forward scattering) to have a charge density wave with both spin and charge gaps are 
consistent with the critical coupling strengths for dimerization obtained analytically and by the quantum Monte Carlo 



method. 14 Away from half filling, the lattice is not statically distorted for purely one-dimensional phonons with 
finite frequency. More realistically, three-dimensionality of the crystal lattice would effectively freeze a periodic lattice 
distortion, if any. In this paper, we do not consider such a case but study effects of dynamical phonons on the spin 
excitation gap in correlated metals. 

Far away from half filling where the Umklapp process is ineffective, a phase diagram is divided basically into 
a gapless phase with a dominant density-wave correlation as in the doped repulsive Hubbard model and a singlet 
superconductor phase with a spin gap for large enough electron-phonon interaction. Near half filling, we found a 
phase with a spin gap and a dominant charge-density-wave correlation between the two phases and call it a spin-gap 
charge-density-wave (SG-CDW) phase. This phase does not have a charge gap and is metallic with a finite Drudc 
weight. The wave number of the SG-CDW is incommensurate to the lattice so that it should be metallic due to 
a phason mode. (l5) The appearance of the SG-CDW phase is thus reasonable and related to the similar phase in 
Ref. H where an Umklapp process and a phonon-assisted forward scattering are not considered. In this paper, we 
clarify its relation with the Mott insulator phase and with a superconductor phase by treating both electron-electron 
and electron-phonon interactions at and away from half filling. Furthermore, effects of the combination of these two 
kinds of interactions near half filling are demonstrated to depend strongly on the form of electron-phonon interaction. 
It is reminiscent of the two-dimensional case in that a system with strong on-site repulsion is more sensitive to an 
intersite electron-phonon coupling than an on-site one for a small phonon frequency. 

This paper is organized as follows: Sec. || introduces the bosonized model. Relations of electron-phonon parameters 



to the g-ology parameters are examined. Sec. Ill outlines the derivation of renormalization group equations. Some 
modifications to previous works are described. It is confirmed that the antiadiabatic limit is correctly reproduced. 
Sec. [Ey| shows phase diagrams for different electron-phonon couplings and phonon frequencies. The evolution of 
parameters with changing scales is demonstrated for typical examples to clarify the roles played by different processes. 
Sec. summarizes the present work. Part of the results presented in this paper were reported briefly elsewhere. 

II. MODEL 

We consider a Tomonaga-Luttinger model extended to include the spin and charge degrees of freedom, a backward 
scattering, an Umklapp process, and phonon-assisted back and forward scatterings. The electronic part is a standard 
one. 0J^] Assuming weak interactions, one can linearize the electronic dispersion relation around the Fermi points. 
The kinetic part is written as a term bilinear in boson operators for spin and charge densities, 

°j(p) = 7/| (Pjt(p) - Piiip)) » (2- 1 ) 

Pj(p) = ^(pjl(p)+Pjl(p)) , ( 2 - 2 ) 
respectively, with the Fourier components of the electron density operators, 

Po,siP) = C },k+ P .s C 3,k,s , (2-3) 

fc 

for right (j—l) and left (j=2) movers with spin s. They obey commutation relations, 

M-P), vitf)] = h(p), M-P')\ = f^<W , (" = <r,p). (2.4) 
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A backward scattering between parallel spins (gi\\), a forward scattering which couples the two branches, j = 1,2, 
(.92 = 92\\ = .92_l), and a forward scattering within a single branch (174 = g^±, g 4 || is set at zero.) are also written as 
bilinear terms. By combining them with the kinetic part, the spin and charge parts of the bilinear terms read 



H„,o 



2-kvp - gi 



^WiiP^ii-p) +er2(-p)er 2 (p)] - -jj-^2[a 1 (p)a 2 {-p) + <Ti(-p)<T2(p)] 

p>0 



p>0 



(2.5) 



Hpfi = 



U [Pl(P)Pl(~P) + P2(-P)P2(P)] - — j ^ [Pl(p)Pl(-P) + P1(-P)P2(P)} 



p>0 p>0 



(2.6) 



respectively, where v-p is the Fermi velocity. 

The other g processes are written with the help of the boson representation of the right-going (j = 1) and left-going 
(j = 2) electron fields with spin s (s = ±1), 



ipj s (x) = lim 



±ikFX =F i-^= (4>p{x) + s<j>„{x)) + i-^= (0 p (x) + s6 a (x)) , for j = j * , (2.7) 
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where a is a cut-off parameter, kp is the Fermi wave number, cj> v (x) and 6 v (x), [y = cr,p), are the fields defined as 



(2- 



(2.9) 



The field H v (x) defined as Il„(x) = ^d x 6 v {x) is conjugate to the phase field (j) v {x), 

[(f> v (x), II M (j/)] = i8 v ^5{x - y) . 

With no magnetic field and at half filling, a backward scattering between antiparallel spins (<7i_l) and an Umklapp 
process (53) read 

ff<r,int = Y / dx V'l.sWV'l ,-,(# -«(a;)^2, s (a;) 

= /o' 9lJ ^2 / ^cos(2V20 (T (x)) , 
(27ra)^ J 



(2.10) 



-ff p ,int = .93 / da; ^Vi, T ( X )V4, 1(^)^2,4 (a#2,t (x) + h.c.) 



2ff3 



da:cos(2v / 2(/)p(x)) , 



(27ra) 2 

respectively. The bilinear terms are rewritten with the fields <j> v (x) and 9 v (x) as 



" / s 



u„A„ (^^(x)) 2 + — ^ (<9 x< ^,(a;)) 2 



with velocities, 



u 2 a = (v F 



9±y _ (in) 2 

2tt) \2tt) ' 



<-(-+£)'-( 



2tt 



and correlation exponents, 



A. 



o _ 27Tt> F ~ ff4 + gl|| 

2ttv f - g 4 - gn\ 



(2.11) 
(2.12) 

(2.13a) 
(2.13b) 

(2.14a) 
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Ki 



2ttv f +g 4 + 3i|| - 232 



p 2ttv f +.g 4 || + 23 2 ' 

With a magnetic field h or away from half filling, we add 

f ^J~2 f V^2~ 

h I dx—d x (j) < ,{x) + fj, / dx—d x <f> p (x) , 



(2.14b) 



(2.15) 



where fi is the chemical potential set to be zero at half filling. In order to work with fixed magnetization (m) and 
filling (Sn), where m = ±1 at saturated magnetization, Sn — at half filling, and Sn — 1 at complete filling, Legendre 
transforms are performed. Kjj Defining 



1p a (x) = 4> a {x) + -=!M , 

V2 



(2.16a) 



ip p (x) = 4> p {x) H — =<5n:c 
v 2 



(2.16b) 



and rewriting ip rT (x) and ipp( x ) as 4'cr(x) an d <j> p (x), respectively, we get the electronic part of the model. The spin 
part, 



H„ — Hi 



0,<T 



1_L 



consists of the "noninteracting" part, 



U ^ (&^(^)) 2 + ^ («9,0 CT (^)) 2 + ^^F^ 



and the perturbation due to the backward scattering 

2$1_L 



(27ra) 2 



cos(2v / 2</> C r(a;) — 1-kuix) , 



while the charge part, 



consists of the "noninteracting" part, 



dx 
2^ 



Hp — Hq.o + H3 



UpKp{d x e p {x)) 2 + (d x <j> p (x)) 



K 



Tr 2 Up5n 2 
2K D 



and the perturbation due to the Umklapp process 
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(2vaf 



dx cos(2y/2(j)p(x) — 2irSnx) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



Here, we do not take the filling dependence of the Fermi velocity into account. In the lowest-order perturbation, the 
Hubbard model with the on-site interaction U corresponds to g^u = g\x = 9i = 33 — 34 — U . g] 

The phonon field is decomposed into a part with momenta near zero, 4>o(p), and the rest with momenta near it (in 
units of the inverse lattice spacing), (f>^(x). The former causes a forward scattering and is represented in the momentum 
space. The latter causes a backward scattering and is represented in the coordinate space. This decomposition is 
therefore made possible by the decomposition of the electron field into the right- and left-going ones so that it is 
achieved unambiguously. Their conjugate momenta are denoted by Ho(p) and H^(x), respectively. 



[<po(p),U (q)] = iS Ptq , 



(2.23) 



h n (x),U n (y)] = iS(x - y) 



(2.24) 



Note that ^(x) and H n (x) are real fields. 
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In the phonon part, we take dispersionless phonons for simplicity. Later we will consider electron-phonon couplings 
in the Holstein model and in the SSH model when evaluating correlation exponents. The neglect of the dispersion of 
bare phonons is not essential since it is absent in the former model and the acoustic phonons are decoupled from the 
low-energy electronic spectrum in the continuum limit of the latter model. fl4|| We divide the phonon part into a part 
with momenta near tt and the rest with momenta near zero, 



H ph = \J dx [IQW+Jffiix)] + [I U o(p) | 2 +^ 2 I MP) I 2 ] , 



(2.25) 



where lj^ and u>o are their respective phonon frequencies. 

The Holstein coupling to the dispersionless phonons is represented by 



V- (a K 2 1 



(2.26) 



with the electron density at site i, rii, the lattice displacement its conjugate momentum pi, the coupling strength 
/?, the spring constant K, and the ionic mass M. In the continuum limit, it can be written as H\ + Hi + if p h, where 
Hi denotes a phonon-assisted backward scattering, 



and H2 stands for a phonon-assisted forward scattering, 

V2 



#2 = 72-p E M-P) + Pi{-p)] 0o(p) , 



(2.27) 



(2.28) 



with 7! = 72 = P/VM and = ljq = \/kJM. Note that the operators tpi. s (x), tp2,s(x), and 4> n (x) above concern 
with the slowly varying part, whereas the rapidly oscillating parts give the factors e ±l7rSnx = e ±»(2kp— ttJx^ 
The SSH coupling to the phonons is represented by 



E 



*(<?»+! - Qi) E( C L C *+M + 4+1 



K 1 
s c i,s) + yfe+i ~ *) 2 + T^fPi 



(2.29) 



with Ci tS annihilating an electron with spin s at site i, a s being the coupling strength, and qi, pi, K and M as above. 
In the continuum limit, by neglecting the phonon dispersion, it can be written as H-± + H p h, where denotes a 
phonon-assisted backward scattering, 



ff-i=7-iE / dx ^LW^-'i 



#1 s (x)^ 2 , s (x)e-^ Snx ] <j> n {x) , 



(2.30) 



with 7_i = 4a s /yM and = 2WK/M. The scatterings i2i and H_i differ in the form factor. 

For later convenience, we define parameters X v depending only on K u and dimensionless coupling strengths Y v as 



X v = 2(1 - K- 1 ) , 



(2.31) 



Ya = 5ij7(t™f) ) 
In the lowest order, the parameters X a and X p are given by 

Xa 5i||/(7TUf) , 

Xp ^ (,9i|| - 2g 2 )/(iTV F ) 



(2.32a) 
(2.32b) 

(2.33a) 
(2.33b) 
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Hp with X a 



Y n 



-X n 



Y, 



Y, 



respectively. Hereafter, the Hubbard model is meant for Hu = H a 
U/(ttvf) > 0, the Holstein- Hubbard model for _ffjj— h = Hn + iJ p h + Hi + H2 with u„ = u>o = u> = J K/M and 
7l = 72 = /3/a/M, and the SSH-Hubbard model for ffssH-H = Hs. + H ph + H-i with = ui = 2y/K/M and 
7_i = Aa s /y/M. Numerical evaluations are performed also in the SSH-Hubbard model supplemented with a forward 
scattering, HgSH-H-f.s. = -ffssH-H + H 2 , with lo = and 72 = 7_i. 

Effective interactions are obtained through the integration of phonon fie lds. First we consider the antiadiabatic 
limit (u> — ► 00), where the interactions are instantaneous. In this limit, Eq. ( 2.26| ) becomes 



which corresponds to g x \\ = g\±_ = g2 = ,93 = gi = —/3 2 /K, and Eq. ( 2.29 ) becomes 



2K 



( c L c i+l,s + c i+l,* c <,»)( c i, a ' c i+l 1 «' + 4+1,8'*,*') 



(2.35) 



which corresponds to gi« = gi± = —33 = — Aoq/K and g 2 = 34 = 0. Later we will take either a phonon-assisted 
backward scattering originating from the on-site coupling (71) or that from the intersite coupling (7-1), but not both. 
We then define a dimensionless coupling strength Yi as 



Y x = 



ill 

TTVfOJ 2 



If necessary, the two cases above are distinguished by / 
coupling strength Y 2 is defined as 



1 for 71 and / = — 1 for 7. 



Yo 



(2.36) 

Another dimensionless 
(2.37) 



The subscripts of Y\ and Y 2 stand for back and forward scatterings, respectively, as those in g parameters. In the 
Holstein-Hubbard model and the SSH-Hubbard model supplemented with a forward scattering, we define Y p ^ as 
Yi = Y 2 = Yph. Summarizing the antiadiabatic case, electron-phonon interaction shifts the electronic par ameter s 



by X n -> X „ - Yu Y„ -> Y„ -Y u X p -> X p - {Y 1 - 2Y 2 ), and Y p -> Y p - fY v From the definitions, (2.32a) 



( |2.33a|) , and (|2.33b| ), the first three are interpreted as due to effective attraction. The difference between the on-site 
coupling (/ = 1) and the intersite coupling (/ = — 1) appears in the effect on Y p . Its origin is easily understood in the 
continuum limit. The effective gi processes come from the contraction of ipt s ( a 0V'i,«( a; ) and W s > (x)ip 2 s ' (x) in Hi or 
H-i, whereas the effective g% processes come from the contraction of ip\ ,(i)fi,j(!c) and ip\ _ s (x)ipi ! - s (x) and that 

of ip\ s (x)ip 2 ,s(x) and ip\ _ s (x)il> 2t - s (x), giving the sign factor /. Hereafter, we will not consider a 34 process, which 
only yields a renormalization of the Fermi velocity. 

To treat finite phonon frequencies, we consider the partition function Z = Tie~" H at the inverse temperature (3 
and map the present system to a two-dimensional classical system in the Euclidean space r = (x,y = vft) with an 
imaginary time r. [pj Because the phonon fields are bilinear in the Hamiltonian, they are analytically integrated out. 
Effective interactions Hi <e g (i = ±1,2) are defined as 









f P 


T r exp 


- / dTH iyCS (T) 


= (^- T ex P 


- / drHiir) 




Jo 




Jo 



(2.38) 



P h 



where T r is the time-ordering operator in the imaginary time, and (• • -) p h denotes an expectation value averaged over 
phonon fields. It is straightforward to derive, in the f3 — ► 00 limit, 



#l,eff(T) 



7f 



cj2(27ra) 5 



ei=±l 



dx / dAr 



ZL e -av|Ar| 



- exp jiei (v^M 1 ) 



<=2,<^ ,e'=±l 



ie 2 ^V / 20 CT (2) — nmx^j + ie'i (^/2(f> p (l) — nSnx^j — ie 2 ^V2(f> p (2) — nSnx^j | , 



(2.39) 



h 2 Mt) = -^t 



L ^ 

p>0 • 



/ dAr 


"i£0 u, |Ar|' 




.2 J 



{pi(p,n) + p2(P)n)}{pi(-p,T 2 ) + /?2(-p,T 2 )} 



(2.40) 



G 



fl-i,eff(r) 



7*1 



^(27ra) 2 



E 

ei=±l 



dx / dAr 



r|Ar| 



E 

e 2 ,<,4=±l 



sgn(eie' 2 )exp 



-ie 2 



(V2(f> (T (2) - nmx) + ie[ W2<j> p (l) - n8nx^ - ie' 2 \V2(p p (2) - irSnx^j i , 



(2.41) 



where arguments 1 and 2 denote (a;, vfTi) and (x, i>fT2), respectively, r = (ri +T2) /2, an d At — ti — r 2 . In the uj — > oo 
limit, [• • •] above becomes <5(ti — r 2 ) so that the e\e 2 — e[e 2 = 1 terms in Eqs. (2.39) and Q2.41 ) correspond to the 
gi\\ process, —e\e 2 = e[e 2 = 1 to the gi± process, £162 = — e[e 2 — 1 to the <fc(= ff3_i_) process, and — eie 2 = — e[e 2 = 1 
to the 03|| process not considered here, whereas the cross terms in Eq. ( |2.40 ) correspond to the g 2 (= g 2 \\ = g 2 ±) 
processes, and the rest to the g A \\ = g t ± processes. Note that =F#4 = T94± in Eqs. ( |2.5| ) and (2.6) are replaced by 
.94|| T 94± if 34H 7^ 0. Again, we can derive X a ->■ X a - Yi, Y CT -> Y a - Yi, X p -> X p - (Yi - 2Y 2 ), and Y p ->Y p - /Yi 
in this limit. 

Finally, we have the "noninteracting" part, Hq = Hq^+Hq. p , and perturbations, Hi±, H3, H±i, e g(r), and H 2te g(r). 



III. RENORMALIZATION EQUATIONS 

Renormalization group equations are derived in a way similar to Voit-Schulz || and Giamarchi-Schulz [0. Because 
of the similarity, we will not go into details but will give an outline detailed enough to explain a few modifications. 
We consider the correlation function, 

fl„(r 12 ) = ^T T e 4V5 ^ (ri) e- l ^ (r2) ^ , (3.1) 

with 1*12 = i"i — r 2 and develop in powers of Y a , Y p , Yi, and Y 2 . In the zeroth order, it is given by 

R W { r 12 )=^ P (-^lo g ^±^-) . (3.2) 



In the expansion, we sum only over those configurations which satisfy the spin- and charge-conservation conditions: 
the sum of the coefficients of the </> CT terms in the exponents and that of the <f> p terms are zero. Each term is written 
as a product of a spin part and a charge part the exponents of which are linear in CT and <p p , respectively, so that 
the expectation value is easily obtained. As a result, each term is written as an integral the variables of which can be 
regarded as positions of Burgers vectors in the Euclidean space Q and denoted by rj, (i = a, b, ■ ■ ■). The perturbation 
-H±i,eff( T ) produces two vectors each of which is allowed to take one of four possible orientations ±bi and ±b 2 with 

I/2I/2 1/2 

bi,2 = (K a , ±K P ). The perturbation H\±_ produces a vector b^ or —h a with b CT = (2K a , 0), while H3 produces 

1 /2 

a vector h p or — b p with h p — (0,2K P ). The perturbation H 2:e ff(r) is represented with different operators and 
not interpreted as producing Burgers vectors: it will be treated later. For example, the 0(Y^) term, (v = a,p), is 
regarded as arising from a logarithmic interaction, 



b(r a ) • b(r b )- log , (3.3) 

2 or 



with b(r a ) = — b(rb) = ±b„ in a classical system of Burgers vectors if m — and Sn — 0. With m ^ and Sn 7^ 0, 
the corresponding interaction becomes anisotropic in the Euclidean space. [ fl7| The O(Yi) term corresponds to a 
contribution from a logarithmic plus linear potential, 

1 / , , , si, |r a -r b | 2 W \y a -Vb\ ,„ .s 

-b(r a ) b(r b )-log — 2 , (3.4) 

& OL q-n 

with a constraint x a — Xb, b(r a ) = — b(r b ) = ±bi, (i = 1,2), and £„■ = vf/uj^. The second term and the constraint 



come from the phonon propagator in Eqs. (2.39) and (2.41) 



Ar(r« - rb) = ^~S(x a - x b ) cM- V \ ~ ] , (3-5) 

which becomes isotropic in the antiadiabatic limit. The integral in each term is dominated by configurations with r a 
very close to r b . Then it is convenient to change the integration variables to R = (r a + r b )/2 and r = r a r b and 
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expand the integrand to the second order in r. Integration can be performed over R. []17j The remaining integration 
over r is cut off at a (r > a) and the function log[(r 2 + a 2 )/a 2 ]/2 is replaced by log(r/a). 
We will search a renormalized correlation function of the form, 

\ 2 a 2 r( 2 J 

regarding the corrections as modifying the exponent K v to an effective one K„ and producing the anisotropy pa- 
rameter dl s . Changing the cut off a from a(l) = ae l to a(l + dl), a pair of Burgers vectors with a(l) < r < a(l + dl) 
are no more visible. This portion of the integrals can be absorbed into changing K u (l)[— 1/(1 — X v (l)/2)] and d v {V) 
into K v (l + dl) and d u (l + dl). Upon rescaling the remaining integrals so that they again range from a to infinity, 
the effective parameters K^ s and d° s are written with modified couplings Y V {1 + dl) and Y\(l + dl). In addition, 
nonneutral pairs of Burgers vectors with a(l) < r < a(l + dl) appearing in the 0(Y u Yx), 0(Y 2 ), 0{Y 2 Y\) terms, 
etc., are combined into new vectors: t>i and b2 into b^, bi and — b2 into b p , h a and — b2(— bi) into bi(b2), b p 
and b2(— bi) into bi(— b^), etc. They also contribute to modifying Y v (l) and Y\(l) into Y v (l + dl) and Yi(l + dl). 
The magnetic field h and the chemical potential /i are evaluated perturbationally. From the scaling h(l) = e l h and 
/j,(l) = e'/i, [[l7| we der ive the relation bet ween m(l) and m(l + dl) and that be tween 6n(l) and 5n(l + dl). The factors 
e ±i2« mx in Eq e ± t 2K8nx m Eq an d ( W7r /2)e~ w "l AT l in Eqs. d|39[ ) and ( |2~4l| ) make the correlation 



functions anisotropic in the Euclidean space. A resultant change of the y scale relative to the x scale in R v (yxi) can 
be interpreted as a change of the velocity in the v channel: d\d v {l)\jdl = — [K v (l) /2]dlog[u v (l)]/dl in the lowest order. 
Finally, we obtain the renormalization group equations, 

dX a (l)/dl - -Y 2 (l)J (2TTm(l)a) - Y^D^l) , (3.7) 
dY a (l)/dl = -X a (l)Y a (l) - rx(0Ar(0 , (3.8) 
d(2Trm{l)a)/dl = 2Trm{l)a + Y*(l) Ji(27rm(Z)a) , (3.9) 

d\og[u a {l)]/dl = ~K a {l) [Y 2 (l)J 2 (27Tm(l)a) + Fi(0Ar(0] » (3-10) 

dX p {l)/dl = -Y 2 (l)J (27rSn(l)a) - Yi(J)Ar(0 + 2^(0^(0 , (3.11) 
dY p (l)/dl = -X p (l)Y p (l) - fY x {l)D„(l) , (3.12) 
d(2i:dn(l)a)/dl = 2nSn{l)a + Y 2 (l) J 1 (2 7 r5n(Z)a) , (3.13) 

dlog[u p (l)]/dl - -\K p (l) [Y 2 (l)J 2 (2ir6n(l)a) + Y^D^l)] , (3.14) 

dY x {l)/dl - [2 - K a (l) - K p (l) - Y a (l)J (2nm(l)a) - fY p (l)J (2ir5n(l)a)] Y x (l) , (3.15) 

dY 2 (l)/dl = 0, (3.16) 

where D P {1) — [ctQ,)/^] exp[— a(l)/£, p ], (fi = 0, ir). Since the cut-off parameter a should be physically of the order of 
fcp 1 , D^(l) can be regarded as D^l) = [u^/ E(l)] exp[— cj^/EQ)], with E(l) = Epe~ l and Ep being the Fermi energy. 
The J n (x) are the Bessel functions, whose oscillations spuriously come from the use of a sharp cut off in real space. 
In numerical integrations, we replace them by J n (x)0(jo,i— I % I) with tqi being the fi rst zero of Jq(x) [Jo(jo,i) = 0], 
6[x) — 1 for x > and 9(x) = for x < 0. The first terms in Eqs. (3^8) and ( 3.1 2| ) are linearized in X v since we 



deal with the lowest-order perturbations. In the absence of Yi and Y2, these equations are reduced to known ones. Q| 
Furthermore, they are consistent with the isotropy in spin space [X a (l) — Y a (l) when m = 0]. The velocities u v would 



be modified also through Eqs. (2.13a) and (2.131), but it is of higher order and neglected here. In the lowest order 
considered here, they are decoupled from the rest. Compared with the previous works p|jl§, the factor e l is absorbed 
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into the prefactor of D p (l) so that the right-hand sides of Eqs. (3.15) and ( 3.1 6| ) are not (3 — • • -)Y\{1) an( A ^2(0 an Y 
more, respectively. If the function (T r [</v(i"i) — 4> v {y2)] 2 ) were renormalized as in Ref. |J, a factor of 2 would appear 
in th e abo ve equations. To keep the equations similar to those in the references, we doubled the definition of Y\. In 
Eq. ( 3.11 ), the effect of a phonon- assisted forward scattering is twice as large as that of a phonon-assisted backward 
scattering, as is expected from Eq. ( 2.33b| ), but it is not in Ref. [lq|. As we show below, the present equations correctly 
reproduce the antiadiabatic limit discussed before. 



To reproduce the antiadiabatic limit, we should recall that log[(r 



2 )/a 2 ]/2 was replaced above by \og(r/a) before 



the renormalization procedure. Because the largest electron-phonon contribution comes from r ~ (fj, = 0, tt), it 
would be valid if a were much smaller than However, a large value of ui^/Ep corresponds to a large value of 
a/^fj,, so that this replacement is no longer valid. Instead, we will deriv e t he Y \ and Y 2 ter ms in the equations for X v 
and Y u , (y = a,p), in a more straightforward manner. Equations ( 2.39| ), ( [2.4QD , and (2.41) for the effective couplings 
Hi je s, {i = ±1, 2), represent the electron-phonon contribution to tl 
(w M /2) exp(— w M | At |). As the time scale changes from E(l)~ 



re effective g parameters with a retardation factor 
to E(l + dl)^ 1 , the integration of the above factor 



over the time slice with width E(l) 1 dl can be absorbed into definitions of the g parameters, for example, as 



9i\\(l + dl)= SiuG)- 2 ^ 



-E(iy 



-E(l+dl)- 



dAi 



E(l+dl) 



E(iy 



g n (l)-l^D„(l)dl 



(3.17a) 



and similarly as 



9i±(l + dl) = g 1± (l)-^D n (l)dl, 



(3.17b) 



gs{l + dl) = g s {l)-f^D n (l)dl 



(3.17c) 



g 2 {l + dl) = g 2 (l) 



7l 



D (l)dl 



(3.17d) 



and an equation for g±(l + dl) similar to that for g 2 (l + dl). Note that here we do not take account of the effective 
g processes of the purely electronic origin. From these equations, the Y\ terms are derived again and the Y 2 term 
is derived here. Since H 2-c s(t) does not produce Burgers vectors, Y 2 is independent of the scale. It simply modifies 
X p (l) without interference with othe r qua ntitie s. No te that similar renormalization group equations can be derived 
with Feynman diagrams. Equations ( 3.1 5| ) and ( 3.16| ) are interpreted as originating from corrections to the electron- 
phonon vertices. The absence of a correction to the Y 2 {1) term is due to the fact that the bubble diagram in the zero- 
sound channel appearing in corrections to the phonon-assisted backward scattering gives the logarithmic correction 
[Eq. ( 3.1 5| )] while that in the third channel appearing in corrections to the phonon-assisted forward scattering does 
not give a logarithmic correction [Eq. ( 3. 16 )1. 

Returning to the discuss ion of the a ntiadi abatic case , we should note that the variable I ranges from to oo so that 

), ( 3.11 ), and ( 3.12| ). It corresponds to E F 



E(l) from E F to in Eqs. (^7|), (|3. 
lattice dynamics at <| At |< oo. 



yp 1 <| At |< oo in spite of unconstrained 
This is in co ntras t to th e pur ely electronic processes contributing to the effective 
g pro cesse s [the first terms in Eqs. ( |3.7D , (3.8), (3.11), and (3.12)] and th e cor rections to the electron-phonon vertex 
[Eq. ( 3.151 )1, which are limited to Ep > E (l) > [no t included in Eqs. ( 3.17 )1- Therefore, we integrate the second 



and third terms in Eqs. ([Tt]), (|3~8|) , ( |3.1iP , and ( |3.12j ) over -oo < / < corresponding to <| At |< E f 1 with fixed 



values of Yi and Y 2 before integrating the renormalization group equations over < I < oo. It leads us to choose the 
initial conditions as |18 



X a (0)=X a -Y 1 (l- e -^/ £F ) 



(3.18a) 



Y a (0) = Y a -Y 1 (l-e-^' E * 



(3.18b) 



X p (0) =X p -Y 1 (l 



+ 2Y 2 (l - e ~ Uo/EF ^j 



(3.18c) 
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Y p (0)=Y p - fY 1 (l-e-^' E ^ . (3.18d) 

In the antiadiabatic limit, the Y\ and Y% terms are dropped from the renormalization group equations so that electron- 
phonon interaction simply shifts the initial conditions as discussed in Sec. [n| 

IV. RESULTS 

First we consider analytic properties of the equations. The phonon propagator D p (l) has a maximum at E(l) = to p 
and decreases exponentially for E(l) < oj^. This smoothly cuts off the contribution from electron-phonon interaction 
at energy of the order of phonon frequency. Since Y\ also changes exponentially the actual cut off is different from 
u)n. Except for unphysically strong, intersite electron-phonon coupling (in view of the fact that a s te+i — Qi) must 
be smaller than the bare transfer integral in the SSH model), Yi(l)D„(t) finally decreases exponentially at large / as 
well as ^(O A)(0' Then the Y\ and Y 2 terms become negligible so that the equations for the spin part and for the 
charge part are decoupled to be in the familiar forms: effectively electron-phonon coupling is no more retarded and 
the system becomes equivalent to a purely electronic system with modified coupling strengths. 

It should be noted that a spin-charge coupling at finite energy is a general phenomenon regardless of the presence 
of electron-phonon interaction: spin and charge are separated (in general only) asymptotically, i.e., in the limit of low 
energy and long wave length. Some purely electronic interaction may induce a filling dependence of the asymptotic 
spin correlation through a spin-charge coupling and an Umklapp process at finite energy. However, the Hubbard 
model is not the case since the fixed point of K a , K*, is known to be unity for U > regardless of filling. Recall 
that, if both of the spin part and the charge part belong to the Tomonaga-Luttinger class, the charge-density-wave 
correlation is proportional to r~ K °~ K p, the spin-density- wave correlation to r~ K ° ~ K p, the singlet-superconductor 
correlation to r~ K °~ K p , and the triplet-superconductor correlation to r~ K <> ~ K p , in the asymptotic limit. 

Deviation from half filling and nonzero magnetization, if any, grow exponentially with I. If a charge density wave is 
formed at half filling, charged solitons would be formed upon doping to connect two degenerate phases. The quantity 
Sna can be regarded as the expected number of charged solitons inside the length a. As the renormalization proceeds 
with increasing length scale and decreasing energy scale, the Umklapp process Y p becomes ineffective at finite energy 
which increases with doping. If magnetization exists, a similar argument holds: the backward scattering Y a becomes 



ineffective at finite energy which increases with magnetization. Then both terms on the right-hand side of Eq. (3.7) 
are finally integrated out to make X a (l) converge to some value: K a (l) does not vanish so that the spin excitation 
spectrum is gapless. 

Hereafter we always assume that magnetization is zero. Then, for spin-independent coupling, X a — Y a , as in the 
Hubbard model, the two quantities are equal at any energy scale, X a (l) = Y a (l). After electron-phonon interaction 
is integrated out (I > lo), it becomes of the form, X a (l) = [I — Iq + X^ 1 (Iq)] . There are two possible fixed points: 
X* = and K* = 1 if X a (lo) is positive; X* — —00 and K* = if X a (lo) is negative. In the former, the asymptotic 
spin correlation is the same as in the noninteracting case (gi\\ =0). In the latter, the spin part corresponds to the 
Luther-Emery state |l3f so that the spin excitation spectrum has a gap. For u> — 00, the spin gap opens if Yph > Y c \. 
For Y c \ = 0, the spin gap opens whenever Yph > regardless of the phonon frequency since X a (0) is negative and 
X a (l) decreases at an increasing rate. 

At half filling, the Umklapp process remains effective unless it is renormalizcd zero. After electron-phonon inter- 
action is integrated out (I > lo), X p (l) and Y p (l) flow on a line X 2 (l) — Y p 2 {l) = X 2 (l ) — Y p 2 (l ) in the direction of 

decreasing X p (l). If X p (l ) >| Y p {l ) |, the fixed point has X* = JxJ^o) - Y p 2 (Z ) and Y* = so that K* > 1. 
The superconductor correlation is dominant if K* > 1. In the Holstein-Hubbard model with Y p h > l^i and u> = 00, 
the initial conditions X a (0) = 3^(0) = -X p (0) = Y p (0) = F el - Y ph lead to K* = and K * = 1 so that it has 
a spin gap and degenerate singlet-superconductor and charge-density-wave correlations. It is consistent with known 
results for the equivalent system, namely, the attractive Hubbard model. If X p (Iq) <| Y p (Iq) \ on the other hand, 
X P {1) flows to —00 so that K* = and the charge excitation spectrum has a gap. This happens quite generally. For 
u) = 00, it applies to the SSH- Hubbard model regardless of supplementation with a forward scattering and to the 
Holstein-Hubbard model if Y ci > Y ph . For the SSH-Hubbard model (/ = -1 and Y 2 = 0), the relation X p (l) = -Y p (l) 
holds at any energy scale: the charge gap always opens regardless of the phonon frequency since X p (0) is negative 
and X p (l) decreases at an increasing rate. For the Holstein-Hubbard model (/ = 1 and Y 2 = Yi), the charge gap 
appears to open as is indeed numerically supported whenever the phonon frequency is finite: the initial condition is 
on the line X p (0) = —Y p (0) and the electron-phonon interaction deviates the flow to the X p (0) + Y p (0) < region 
even if X p (0) > > Y p (0) because Y\(l) initially increases and overcomes the effect of Y 2 {1) then. 

Consequently, three phases are possible at half filing: a charge-density-wave phase with both spin and charge gaps 
(SCG-CDW) with K* = and K* = 0; a charge-gap phase with a dominant density-wave correlation (CG-DW) 
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with K* = 1 and K* — 0; and a spin-gap phase with degenerate singlet-superconductor and charge-density- wave 
correlations (SG) with K* — and K* = 1. The above facts of the spin and charge correlations at half filling are 
consistent with analytical and quantum Monte Carlo results for Y c \ = 0: O] the SCG-CDW phase always appears 
except for the Holstein model in the antiadiabatic limit where the SG phase appears instead. The SCG-CDW phase 
in the Holstein model has an ordinary charge density wave, while that in the SSH model has a bond order wave, where 
bond lengths and bond charge density are alternated but the charge density is uniform. Both cases have a similarity 
to a band insulator in the sense that they have both the spin and charge gaps. On the other hand, in the CG-DW 
phase, only the charge gap opens. It is similar to the repulsive Hubbard model. Without a logarithmic correction to 
the correlation functions, it is not clear whether the charge-density-wave order is dominant or the spin-density-wave 
order is dominant. 

Phase diagrams at half filling are shown in Fig. [l]. With Y c \ > 0, F p h must be larger than a critical coupling strength 
to have the SCG-CDW phase. The effect of electron-phonon interaction on a spin gap depends sensitively on whether 
the coupling is on-site or intersite, especially when the phonon frequency is low: the critical coupling strength is 
much lower for the intersite coupling than for the on-site coupling. This is due to the sign factor / in the corr ection 
to the electron-phonon vertex, Eq. ( 3.15 ), and in the correction to the effective Umklapp process, Eq. ( 3.12 ). For 
the intersite coupling (/ = —1), the increasing Umklapp process [Yp(/)] strengthens the phonon-assisted backward 
scattering [Yx(0L contrary to the on-site coupling (/ = 1). The Umklapp process is further strengthened by the 
phonon-assisted backward scattering: thes e tw o processes are constructively interfered with each ot her. The growing 
Yx(l) helps the opening of a spin gap [Eq. ( |3.7|) 1 and the growth of the density- wave correlation [Eq. ( |3.1lj )]. As shown 
later, this difference is evident near half filling also: the intersite coupling favors the spin gap and the density-wave 
correlation. It is reasonable in view of the fact that, in the limit of strong on-site repulsion and small phonon frequency, 
the intersite coupling enables alternation of the superexchange and produces a spin-Peierls state no matter how the 
coupling is small. Meanwhile, the on-site coupling must be large enough to compete with the on-site repulsion to affect 
the system. As the phonon frequency increases, however, the difference between these two couplings becomes small. 
In the antiadiabatic limit, the spin gap opens for F p h > Y c \ in both cases. As a function of the phonon frequency oj, 
the critical coupling strength increases with lj for the intersite coupling, while it decreases for the one-site coupling. 

Away from half filling, the fixed point K* is always positive and the charge excitation spectrum is gapless since 
the Umklapp process is integrated out at finite energy. Therefore, four phases are possible: a gapless phase with a 
dominant superconductor correlation (SC) with K* = 1 and K* > 1; a gapless phase with a dominant density- wave 
correlation (DW) with K* = 1 and K* < 1; a singlet superconductor phase with a spin gap (SG-SC) with K* = 
and K* > 1; and a phase with a spin gap and a dominant charge-density- wave correlation (SG-CDW) with K* = 
and K* < 1. The DW and SG-CDW phases have a finite Drude weight, 2K*u* p , so that they are metallic. The DW 
phase appears in the doped repulsive Hubbard model. The SG-SC phase is regarded as a conventional superconductor 
phase since it is accompanied by a spin gap with the help of electron-phonon interaction, although a true long-range 
order is absent. 

Far away from half filling where the Umklapp process is ineffective, a phase diagram is divided basically into the 
DW phase for Y^^Yp^ and the SG-SC phase for Y p h^,Y e i (Figs. || and |J). For finite phonon frequency, the SC 
phase occupies a narrow region between the two [not shown in Figs. ^| (a) and ^| (c)] and may be an artifact of the 
lowest-order renormalization group approach. However, the spin correlation and the charge correlation are separated 
in the asymptotic limit so that there is no reason to exclude the SC phase. A phonon-assisted backward scattering 
tends to open a spin gap, while it enhances the density-wave correlation. To make the superconductor correlation 
dominant, a phonon-assisted forward scattering must overcome the backward scattering. In the continuum limit of 
the SSH-Hubbard model where a coupling with acoustic phonons is absent, the SG-SC phase is made possible by the 
supplementation with a forward scattering. |l8[ ] 

Phase diagrams for the Holstein-Hubbard model are shown in Fig. |^. For very small phonon frequency u> and near 
half filling, the SG-CDW phase appears between the DW and SG-SC phases. In Fig. || (a), as the filling approaches 
the half (5n — 0), the boundary between the SG-CDW and DW phases goes upwards abruptly in a very close vicinity 
of the half filling, although it occurs too close to distinguish the boundary from the ordinate. With Y p h values shown, 
the spin gap does not open at half filling in Figs. |^ (a) and (b). Interestingly, a spin gap is induced by doping in 
these two cases. The phonon-assisted backward scattering [Yi(Z)] is destructively interfered with the Umklapp process 
[Yp(Z)] so that it is not strong enough to open a spin gap at half filling for the parameters used in the figures. Upon 
doping, the destructive interference becomes weak and the spin gap opens. Upon further doping, the effects of Y p (l) 
and Yi(l) are overcome by the phonon-assisted forward scattering [Y2O] m Eq. (3.11) so that X* > 0, K* > 1, 
and therefore the superconductor correlation becomes dominant. As the phonon frequency lj increases, the SG-CDW 
phase disappears from the higher doping side and finally goes away for a phonon frequency at a fraction of the Fermi 
energy. The fragility of the SG-CDW phase in the Holstein-Hubbard model is due to the fact that the cut off of the 
Umklapp process acts both for opening a spin gap and for strengthening the superconductor correlation. Because of 
retardation, the former effect becomes largest at energy of order w, while the latter effect is independent of retardation. 
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This is why the SG-CDW phase is possible only for small uj. The filling dependence becomes weak with increasing w. 

The above arguments become clear when one observes the evolution of renormalized quantities with decreasing 
energy scale shown in Fig. |^ for the Holstein-Hubbard model. The first zero of Jq(x) is about 2.4 so that the Umklapp 
process is cut off at about I ~ 2 [Fig. [| (c)]. Before reaching the cut-off energy scale, X p (l) rapidly decreases [Fig. ^ 
(b)] and Y\(l) also decreases [Fig. ^ (d)] due to the destructive interference with Y p (l), which is not shown but grows 
exponentially. Since the phonon frequency is very small here, the electron-phonon interaction becomes effective at a 
much smaller energy scale. The products of an electron-phonon coupling and the corresponding phonon propagator 
have a peak around I ~ 7 and decrease exponentially after the peak [Figs. [I| (e) and (f)]. Before reaching I ~ 7, 
X a {l) decreases proportionally to the inverse of I but is still positive [Fig. ^ (a)]. The magnitude of the slope of X a (l) 
becomes large again around I ~ 7 and finally X a {l) becomes negative: it flows to — oo. Because of the destructive 
interference, Y\{1) is smaller than Y\ around I ~ 7, and consequently Yi{l)D T {l) is smaller than Y 2 (l)D (l). Therefore, 
X p (l) increases after the Umklapp process is cut off. However, the electron-phonon effect is not strong enough to 
make X p {l) positive. 

Phase diagrams for the SSH-Hubbard model supplemented with a forward scattering are shown in Fig. ^. Now 
Y p (l) and Y\(l) are constructively interfered with each other. The tendency for opening a spin gap near half filling is 
much stronger than in the Holstein-Hubbard model, regardless of supplementation with a forward scattering. Note 
the difference between the scales of Figs. || and ||. Contrary to the Holstein-Hubbard model, as the filling deviates 
from the half, the constructive interference becomes weak as well as the tendency for opening a spin gap. The 
critical coupling strength for a spin gap is therefore an increasing function of | 5n | , whose slope becomes gentle with 
increasing u). Since Y p (l) and Yi(l) are enhanced by each other, larger deviation from half filling is necessary for 
YqQ) to overcome these two to make the superconductor correlation dominant, compared with the Holstein-Hubbard 
model. As a consequence, the SG-CDW phase appears in a much wider region. As the phonon frequency uj increases, 
the interference effect becomes weak and larger electron-phonon coupling is necessary for a spin gap. Therefore, the 
SG-CDW phase is shifted to larger Y p h values with increasing uj. However, this phase remains in the uj — > oo limit. 
In this limit, the condition for a spin gap is the same as in the Holstein-Hubbard model, while the condition for a 
dominant superconductor correlation is different due to the sign factor / in the initial condition for Y p (0). 

The evolution of renormalized quantities with decreasing energy scale is shown in Fig. | for the SSH-Hubbard 
model supplemented with a forward scattering. Before the Umklapp process is cut off at about I ~ 2 [Fig. || (c)], 
X p (l) rapidly decreases [Fig. || (b)]. Note that F p h is set at 0.1 in Fig. hi while it was set at a much larger value 
of 0.4 in Fig. i. Because of the constructive interference with Y p (l), Yi(l) increases when the Umklapp process is 
effective [Fig. 5 (d)]. Before the electron-phonon effect becomes largest at about I ~ 5 [Figs. || (e) and (f)], X a (l) 
decreases roughly proportionally to the inverse of / [Fig. [s] (a)]. The magnitude of the slope of X a (l) becomes large 
again around I ~ 5 and finally X a (l) becomes negative and flows to — oo. Compared with Fig. ^, the products of an 
electron-phonon coupling and the corresponding phonon propagator grow more rapidly and are not negligible when 
the Umklapp process is cut off. Therefore, the region where the electronic backward scattering decreases X a (l) and 
the region where the phonon-assisted backward scattering decreases X a (l) are not clearly distinguishable. Because of 
the constructive interference, Yi(l) is about twice as large as Y\ around I ~ 5, and consequently Yi(l)D n (l) is about 
twice as large as Y2 (l)Dg(l). Therefore, the effects of the back and forward scatterings assisted by phonons cancel 
each other [Eq. 3.1l| so that X p (l) does not change so much after the Umklapp process is cut off. 



V. SUMMARY 



We have studied asymptotic correlation functions of ncarly-half-filled one-dimensional conductors coupled with 
dispersionless phonons using the boson representation of the electron fields and the lowest-order renormalization group 
method. On-site repulsion is considered for the electronic scattering processes. As the renormalization proceeds, 
the Umklapp process is initially effective even away from half filling and this causes the filling dependence of the 
correlation functions. We consider a phonon-assisted forward scattering which enables the superconductor correlation 
to be dominant. In addition, we take two types of phonon-assisted backward scatterings originating from the Holstcin 
model and from the SSH model. They couple with both spin and charge until the energy scale reaches the phonon 
frequency: variation of filling modifies the spin correlation also. The variety of phases comes from the cut-off energy 
scale of the Umklapp process relative to that of electron-phonon interaction, which depends upon doping. The 
variety also comes from the interference of the Umklapp process and the phonon-assisted backward scattering, which 
is destructive for the Holstein-type on-site electron-phonon coupling and constructive for the SSH-type intersite 
electron-phonon coupling. 

At half filling, a charge gap opens by the Umklapp process in most cases. A spin gap also opens if the phonon-assisted 
backward scattering is strong enough. Since these two processes interfere constructively for the SSH coupling, the 
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critical coupling strength for a charge density wave with both spin and charge gaps is weaker for the SSH-Hubbard 
model than for the Holstcin-Hubbard model. The difference between these two models is large for small phonon 
frequency since the interference is effective for a wider energy range of the renormalization process. This interference 
effect appears also near half filling, where the Umklapp process is integrated out at finite energy so that the charge 
excitation spectrum is always gapless. The spin gap opens more easily by the intcrsite coupling than by the on- 
site coupling, especially for small phonon frequency. The above two processes enhance the density-wave correlation. 
To make the superconductor correlation dominant, the phonon- assisted forward scattering has to overcome them. 
It is achieved upon doping much easier for the Holstein-Hubbard model since the two negative processes interfere 
destructively 

Far away from half filling, a phase diagram is divided basically into the gapless metallic phase with a dominant 
density-wave correlation, if the on-site repulsion is stronger, and the conventional superconductor phase with a spin 
gap, if the electron-phonon interaction is stronger. Near half filling, we found a spin-gap metallic phase with a 
dominant charge-density-wave correlation between the two phases. This phase appears in a wider range for the SSH- 
Hubbard model supplemented with a forward scattering because of the constructive interference as mentioned above. 
For the Holstein-Hubbard model, it appears only for small phonon frequency, where the spin gap is induced by doping 
since the destructive interference is weakened as the filling deviates from the half. 

For the spin-gap metallic phase, one-dimensionality is essential since the spin excitation and the charge excitation 
are decoupled in the asymptotic limit. The gi± process is effective only for the spin excitation and produces a gap if 
strong. For higher dimensions, a considerably good nesting would be necessary for a similar phase. 
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FIG. 1. Phase diagrams of (a) the Holstein-Hubbard model and (b) the SSH-Hubbard model supplemented with a forward 
scattering, at half filling Sn = and without magnetization m = 0. The electronic parameter is set at Y e \ — 0.4. 

FIG. 2. Phase diagrams of the Holstein-Hubbard model for (a) cu = 10~ 3 _E F , (b) cu = 10~ 2 E F , (c) uj = 10 _1 £ F , and (d) 
lj = oo, for Sn / (See Fig. 1 for Sn = 0.) and m = 0. The electronic parameter is set at Y e \ = 0.4. 

FIG. 3. Phase diagrams of the SSH-Hubbard model supplemented with a forward scattering for (a) uj = 10~ 2 -Ef, (b) 
lj = 10 _1 -Bf, (c) u) = E F , and (d) ui — oo, for Sn / (See Fig. 1 for <5n = 0.) and m = 0. The electronic parameter is set at 
Yd = 0.4. 



FIG. 4. (a) X CT (/), (b) X p (l), (c) 27r5n(/)Q, (d) Yi(/), (e) Y 1 (l)D„(l), and (f) F 2 (/)£>o(0, in the SG-CDW phase of the 
Holstein-Hubbard model. Parameters are u> = W~ 3 E F , Y pb = 0.4, 2-irSna = 0.1, m = 0, and Y c i = 0.4. 



FIG. 5. (a) X a {V), (b) X p (l), (c) 27r5n(/)Q, (d) Yi(Z), (e) Fi(/)^(0, and (f) F 2 (0A)(0> in the SG-CDW phase of the 
SSH-Hubbard model supplemented with a forward scattering. Parameters are u> = 10~ 2 Ef, Y p h = 0.1, 2irSna = 0.25, m = 0, 
and Yd = 0.4. 
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